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ABSTRACT 

We investigate the effects of gravitational lensing in the binary pulsar system J0737-3039. Current 
measurement of the orbital inclination allows the millisecond pulsar (A) to pass very close (at i? m i n rs 
4000 km) in projection to the companion pulsar (B), with i? m i n comparable to the Einstein radius 
(rs 2600 km). For this separation at the conjunction, lensing causes small (about 10%) magnification 
of the pulsar A signal on a timescale of several seconds, and displaces the pulsar image on the sky 
plane by about 1200 km. More importantly, lensing introduces a correction (of several fj,s) to the 
conventional Shapiro delay formula used in pulsar timing analysis, and gives rise to a geometric 
time delay together with the delays associated with the pulsar spin period. These lensing effects 
can influence the determination of the system parameters by both timing and scintillation studies. 
Given the current uncertainty in the orbital inclination, more extreme manifestations of lensing (e.g. 
magnification by a factor of several) are possible. We compare our predictions with the existing 
observations and discuss the possibility of detecting gravitational lensing signatures in the system. The 
anomalously high point in A's lightcurve close to superior conjunction might be caused by gravitational 
lensing. 

Subject headings: pulsars: general — stars: neutron — pulsars: individual (J0737-3039A, J0737- 
3039B) — gravitational lensing — binaries: general 



1. INTRODUCTION 

The binary pulsar system J0737-3039 consists of a mil- 
lisecond pulsar (pulsar A with a period Pa = 23 ms) and 
a normal pulsar (pulsar B with Pb — 2.8 s) in a 2.4 hr or- 
bit with eccentricity e = 0.088 (Burgay et al. 2003). The 
system has a nearly edge-on orientation with respect to 
our line of sight, leading to a brief (about 27 s) eclipse of 
pulsar A by the magnetosphere of pulsar B during each 
orbit at A's superior conjunction (Lyne et al. 2004; Kaspi 
et al. 2004; McLaughlin et al. 2004). The depth' of the 
eclipse is quite significant, corresponding to an absorp- 
tion optical depth of order a few (Kaspi et al. 2004). The 
eclipse is asymmetric: it is deeper and longer after the 
conjunction, indicating that the optical depth increases 
from the eclipse ingress to egress. 

Timing observation (including measurement of Shapiro 
delay) provided the first determination of the system's 
inclination angle with respect to our line of sight, i = 
87° ±3° (Lyne et al. 2004). Recently, Coles et al. (2004) 
have used the correlation of the interstellar scintillations 
of both pulsars to obtain a remarkable constraint on 
the minimum distance between A and B on the plane 
of the sky, R m i n = 4000 ± 2000 km, corresponding to 
i = 90.29° ± 0.14°. As we show in this letter, such a 
small projected separation of A and B gives rise to inter- 
esting gravitational lensing effects during each orbit. 

2. GRAVITATIONAL LENSING IN J0737 

At the superior conjunction of pulsar A, the distance 
between A and B projected to our line of sight is an = 
a\ sin i | ( 1 — e 2 )/(l + esinw), where a — 8.79 x 10 5 km 
is the semimajor axis of the orbit, lo is the longitude of 
periastron (at present u> = 74°, and Co = 16.9°/yr). With 
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pulsar B mass Mb = L25 M Q , the Einstein radius is 

R E = {2R g a\\) 1/2 ~ {2R g a) 1/2 = 2550 km, (1) 

where R g = 2GMb/c 2 = 3.69 km. This is already com- 
parable to the minimum projected A-B distance i? m i n = 
a|||tani| _1 ~ 4000 ± 2000 km. Given the present un- 
certainty of the measurement, the real i? m i n might be 
smaller than Re, increasing the strength of gravitational 
lensing. 

The typical timescale of lensing is io ~ -Re/v Rs 4 s 
(where v rs 680 km s _1 is the relative velocity between 
the two pulsars), considerably shorter than the total du- 
ration (~ 27 s) of the magnetospheric eclipse. Gravita- 
tional lensing causes variation of the pulse flux during the 
eclipse, with the magnification factor (Paczynski 1996) 



u\Ju l + 4 -Ke 

Here R s — r(l — sin 2 i sin 2 t/;) 1 / 2 is the separation be- 
tween the two pulsars in the plane of the sky, with 
r = a(l — e 2 )/(l + ecoscp) and i\> = <j) + lo being the 
physical separation between A and B and the longitude 
from ascending node, respectively (<fi is a true anomaly 
measured from periastron). Let ip = ir/2 + Aip (the A-B 
conjunction corresponds to Alp = 0) and i = tt/2 + Ai, 

we have R s ~ ay [(Ai) 2 + (A^) 2 ] 1/2 for |Ai|, |A^| < 1. 

For the current best estimate of i? m in — a»Ai ~ 
4000 km (i = 90.29°) one obtains A roax = 1.13. 
Such 10% flux variation due to lensing could have been 
marginally detectable with currently available observa- 
tions (Kaspi et al. 2004) at 820 MHz if the superior 
conjunction of pulsar A were not accompanied by its si- 
multaneous eclipse by the pulsar B's magnetosphere (see 
W2.2t> . This may however be possible when more data on 
the eclipse become available. 
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Fig. 1. — The combined geometric delay and gravitational 
(Shapiro) delay of pulsar A's signal (for the "+" image) includ- 
ing the lensing effect by pulsar B (solid lines). Timing perturba- 
tions associated with the pulsar spin period (see the text) are not 
included. The dashed lines show the Shapiro delay when lensing 
is not included. Both the solid and dashed lines correspond to 
^min = 0, 1000, 2000 and 4000 km from top to bottom. The dot- 
ted lines show the time delay associated with the "— " images (with 
i? m ; n increasing from bottom to top). Note that the "— " images 
are usually much fainter than the "+" images. The longitude is 
measured from the superior conjunction of pulsar A (i.e., when A 
is behind B). 



Smaller i? m i n leads to stronger magnifications: for 
R min = 2000 km (i = 90.14°), corresponding to la 
deviation from the mean, we find A max = 1.6; for 
i? min = 1000 km (i = 90.07°), A max = 2.8. Thus, the 
amplification of pulsar A signal by a factor of several is 
not inconsistent with the currently measured inclination 
of the system — in the end even the purely edge-on ori- 
entation of the system (i = 90°) deviates only by 2<r from 
the measurement of Coles et al. (2004). 

2.1. Time Delay 

Gravitational lensing has important effect on the delay 
of the pulse arrival time (Schneider 1990). The geometric 
time delay is given by 



(At) 



Rg ( AR± 
c V Re 



(3) 



where AR± = \R± — R s \ is the relative distance in the 
plane of the sky between the image position of the source, 
R± = 0.5(R S ± \JR% + ^R 2 E ), and its fiducial position in 
flat space, R s ; "+(—)" refers to the image which is closer 
(further) from the source. The magnification factor of 
the image is A± = (A ± l)/2. For u = R s /Re -> 
one finds AR± — > Re and the two images merge into an 
Einstein ring; for u 1 one has AR + — ► R^u^ 1 and 
Ai?_ — > RgU, but the "— " image contributes negligible 
light (Paczynski 1996). For J0737, the maximum possible 

geometric delay is (At)^om = R g /c = 12 (is, achieved 
for u = (i = 90° and i> = 90°). 

Gravitational (Shapiro) delay amounts to as much as 
100 (is in J0737. If lensing is neglected, the Shapiro delay 
is given by the standard formula (Blandford & Teukolsky 



(Af)i- l -) = ^lnf 1 l + eC ° S ^ 
v /glav c \l - smismifj 



(4) 



Note that for i = 90°, this expression diverges at tp = 
90° since in the absence of lensing the ray would pass 
through the infinitely deep potential of pulsar B. This is 
an artifact of neglecting the light bending caused by the 
lensing — the correct formula including the lensing effect 



is 



(Ai) grav =-- ^ln 



Rl 
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where r \\ = r sin i sin tp is the separation of A and B along 
the line of sight. It is easy to see that for Re = 
(no lensing), this expression reduces to the usual (no 
lens) result. Unlike eq. the lensing-corrected ex- 
pression never diverges: the maximum possible delay 

is (At)| 1 " a a v x) = (Rg/c) ln[o(l - e 2 )/R g ] = 152 (is, again 
achieved for i = ip = 90°. 

Finally, light bending causes radio pulses detected on 
Earth to be emitted in a direction slightly different from 
our line of sight, by AR+/rn for the "+" image. This 
corresponds to a different spin phase of the pulsar at the 
moment of pulse emission than the phase would be in the 
absence of lensing, leading to another timing perturba- 
tion proportional to the pulsar spin period (Schneider 
1990). Near the superior conjunction of pulsar A this de- 
lay amounts to several (is for the pulsar A signal 5 . This 
type of delay is very interesting since it is sensitive to 
the orientation of the pulsar spin axis (Doroshenko & 
Kopeikin 1995). Its contribution is not symmetric with 
respect to the conjunction point which makes possible its 
separation from other timing contributions. Because of 
the importance of this delay for constraining the system's 
geometry we postpone its detailed study to a separate 
paper (Rafikov & Lai 2005). 

In Figure n we compare the time delay with and with- 
out the geometric delay and corrections to the Shapiro 
delay caused by lensing effects in J0737 for the domi- 
nant, "+" image near the point of eclipse of pulsar A for 
different values of i? m i n - We do not include the delays as- 
sociated with the pulsar spin period in this Figure (this is 
equivalent to setting pulsar spin period to zero). One can 
see that for i? m i n < 4000 km, (Ai) gcom + (At) grav differ 
from the conventional Shapiro delay formula by > 4 (is. 
This emphasizes the importance of accounting for the 
lensing effects which affect the estimate of the system's 
parameters. Both the geometric delay and the correction 
to the Shapiro delay are localized near the conjunction 
of the pulsars and their variation has a typical timescale 
of to (unlike the Shapiro delay without lensing, which 
varies on the scale of the orbital period of the system). 

For R s ~ Re, assuming that pulsar A is a point source 
(see H2.3|l . one should expect to see a splitting of the 
original pulsar A pulse profile in two because of the two 

3 The formula given by Schneider (1990) is incorrect; e.g., it does 
not reduce to the standard formula in the no-lensing limit. 

4 There is also additional contribution to the delay due to the 
rclativistic aberration of light (Smarr & Blandford 1976). 

5 The corresponding delay of the signal of pulsar B at its superior 
conjunction is ~ 10 s (because of the large Pb)- 
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Fig. 2. — Model lightcurve of the eclipse with the effect of 
the gravitational lensing folded in. Modulation of the lightcurve 
by pulsar B's rotation is neglected (see text). The solid curve 
represents unmagnified eclipse profile. Other curves correspond 
to -R m in = 4000 km (short- dashed), 1000 km (dot-dashed) and 
450 km (long- dashed) (we assume that the two images suffer the 
same amount of attenuation; see text). The star approximately 
marks the position of the flux anomaly during one of the eclipses 
reported by Kaspi et al. (2004). 



gravitationally lensed images, so that the modified pulse 
profile is a sum of the two weighted replicas of the nor- 
mal pulsar A pulse profile slightly shifted in time with 
respect to each other. The relative time delay between 
the two contributions and their amplitudes should vary 
in the course of the eclipse as R s varies. This modi- 
fication of the pulse profile may be difficult to detect: 
for i? m i n = 4000 km, one finds that at the conjunction 
At = (At) gcom -|-(A£) grav ~ 138 [is and 218 fis for the two 
images, with the image magnification factor A + = 1.07 
and A^ = 0.07. Thus, the contribution to the total pulse 
profile coming from the "— " image with longer time de- 
lay (see Fig. 1) is very small. Also note that the "— " 
image samples the deeper interior of the pulsar B mag- 
netosphere than the "+" image does, which causes its 
stronger absorption. Nevertheless, the two-image time 
delay corrections described above can in principle be de- 
tected with exquisite timing observations, especially if 
-Rmin is small (e.g., for i? m i n = 1000 km, the time de- 
lays of the two images are At ~ 156 (is and 175 /is, with 
A + = 1.9 and A_ = 0.9; see Rafikov & Lai 2005 for more 
details). 

2.2. Has Lensing Been Detected in J0737? 

Detecting any signature of lensing is made difficult in 
the case of J0737 by the strong flux attenuation of pulsar 
A at its superior conjunction. Magnetospheric eclipse 
affects both the magnification measurements and (to a 
lesser extent) timing measurements. 

Because of the short lensing timescale compared to the 
eclipse duration, the lensing magnification should mani- 
fest itself as a spike in the lightcurve of pulsar A in the 
inner part of the eclipse profile. Using the published 
parameters of the eclipse (Kaspi et al. 2004) we model 
pulsar A lightcurve in Figure at 1400 MHz for differ- 



ent values of i? m i n 6 . The shape of the spike is skewed 
because the magnification profile is convolved with the 
eclipse lightcurve which exhibits an increase in optical 
depth with time. This causes the spike to occur slightly 
before the superior conjunction of pulsar A. 

It is rather curious that strong fluctuations in the flux 
of pulsar A during the eclipse have already been de- 
tected. Kaspi et al. (2004) report a detection of a ha 
outlier in the eclipse lightcurve at 1400 MHz just sec- 
onds before the conjunction (see their Figure 2). The 
same phenomenon is also seen in the measurements by 
McLaughlin et al. (2004) who demonstrate the variation 
of the eclipse profile with the rotation phase of pulsar 
B. In their lightcurves, among the multiple spikes and 
dips caused by B's magnetosphere, an anomalously high 
flux point has again been seen prior to the conjunction 
(see their Figures 1 and 3). Gravitational lensing can 
explain the outlying point in the Kaspi et al. measure- 
ments if i? m in — 400 km (i = 90.029°), corresponding to 
the magnification of ^4 ma x = 6. 

The lensing interpretation of the observed flux anoma- 
lies during the eclipse is complicated by the fact that 
spikes are not observed during all eclipses and there is 
only one strongly outlying point. Also, in the rota- 
tion phase-resolved measurements of McLaughlin et al. 
(2004), one should expect to see a flux enhancement of 
A at the B pulse phase 0.25 similar to that seen at the 
phase 0.75, since in both phases the unmagnified fluxes 
are quite significant. The error bars on all these measure- 
ments are currently too large for us to decide whether the 
flux anomalies are truly caused by gravitational lensing. 

Current timing measurements of pulsar A have not 
been able to detect time delay during the eclipse, only 
the upper limit on the change in the dispersion mea- 
sure during the eclipse has been set at the level of 
D < 0.016 pc cm -3 (Kaspi et al. 2004). This corresponds 
to a current timing resolution of about 10 — 20 fis which 
makes the geometric delay and lensing correction to the 
Shapiro delay hard to observe at present — the effects 
are of order 4 fis for i? m i n = 4000 km (see Fig.^). How- 
ever, with the improved timing resolution at the level of 
1 fis, one would be able to detect them even if i? m i n is as 
large as 10 4 km (at this i? m i n magnification is too small 
to be detectable) . Timing can be a much better probe of 
lensing than the magnification if i? m i n > Re- 

2.3. Size of the Emission Region of Pulsar A 

Under many circumstances lensing can be used as a 
probe of the size of the source (Paczyhski 1996). The size 
of the emitting region in pulsar A is limited from above 
by the light cylinder radius, i? emj A < Rl,a — c/^A ~ 
1000 km. Finite source size limits the magnification to al- 
ways stay below ^/4(i? £ /i? cm!A ) 2 + 1 ~ 5{Rl,a/ 'Rem,A) 
(e.g., Witt & Mao 1994), and any observed deviation 
from the point-source amplification formula could con- 
strain i? em A. But unless R mm < Rem, A < 1000 km this 
limit cannot be achieved and the finite source size effects 
are small. If the emitting region size in the millisecond 

6 Note that the lightcurve of Kaspi et al. was based on aver- 
aging over 2-second time intervals of the pulses. On a finer time 
resolution ( 0.3 s), the lightcurve shows many spikes and dips due 
to the modulation of A's signal by B's rotating magnetosphere 
(McLaughlin et al. 2004). 
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pulsar is much smaller than Rl.a, the detection of the 
source size through lensing would be unlikely. 

2.4. Effects of Orbital Evolution 

Effects of general relativity are very strong in the J0737 
system. They drive the evolution of its orbital param- 
eters and affect its geometry, which directly translates 
into the change of the lensing signal. Geodetic precession 
causes a wobble of the orbital plane of the system on a 
period of 75 yrs (Burgay et al. 2003). Unfortunately, the 
amplitude of the inclination variation caused by this is 
very small, of order (SU/X or b) sm ^S£ ~4x 10~ 5 sin#s£, 
where Sa is the spin angular momentum of pulsar A, 
L or b is the orbital angular momentum, and 9s l is the 
angle between the spin axis and the orbital angular mo- 
mentum axis. Thus, geodetic precession is not going to 
affect gravitational lensing in J0737 unless |Ai| < 0.02°. 

More important would be the effect of the periastron 
advance on a timescale of 2f yrs. It causes the projected 
distance between the two pulsars at their conjunction, 
-Rmin, to vary from a(l — e)| cosi| to a(l + e)\ cosi|. Thus, 
in the course of the apsidal precession, i? m i n varies by 
about 18%. This causes a change in the lensing signa- 
tures (both magnification and timing) at a similar level 
of ~ 10% and might be detectable in the future. 

3. DISCUSSION 

Additional propagation effects are caused by the 
plasma refractivity in the magnetosphere of pulsar B dur- 
ing the eclipse (e.g. Emmering & London 1990). Assum- 
ing that the magnetospheric charge density is / times the 
Goldreich- Julian density (Goldreich & Julian 1969), the 
deviation of the index of refraction from unity is of order 
An ~ 10~ 7 / / v 2 , where v is the radio wave frequency in 
GHz and density is evaluated at the Einstein radius Re 
7 . By contrast, the effective index of refraction due to 
lensing deviates from unity by R 9 /Re — 1-5 x 10~ 3 . For 
the J0737 system, / is likely of order 100 (Rafikov & Gol- 
dreich 2005), in which case the propagation effects due to 
the plasma in the vicinity of pulsar B are not important 
compared to the gravitational lensing effect. 

Since the gravitational lensing signal is strongly depen- 
dent upon the minimum projected separation i? m ; n of the 
pulsars in the plane of the sky, the detection of lensing 
effects could provide tighter constraints on the system's 
inclination, which is currently known at an accuracy of 
0.14°. Timing would be a good probe of the inclination 
if i? m i n > Re while the magnification signal would also 
be useful for handling the case i? m in ^ Re- 

We have already mentioned that the correction to the 
Shapiro delay and the additional geometric delay caused 
by lensing at the binary conjunction may affect the es- 
timation of the system's parameters. This suggests that 
lensing must be accounted for when dealing with the ef- 

7 This estimate assumes cold plasma in the magnetosphere. For 
a rclativistic plasma, An is reduced further, making the refractive 
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fects of general relativity in J0737. Gravitational light 
bending is also relevant for the measurements using the 
correlation of scintillations of the two pulsars. The anal- 
ysis of Coles et al. (2004) assumes that pulsar A's tra- 
jectory in the plane of the sky is a straight line at the 
conjunction. This is not true if lensing is strong because 
it displaces the position of pulsar A on the sky along a 
curved path deviating from the straight line trajectory 
by AR + . Even for i? m i n = 4000 km, the shift can be 
as large as 1200 km, almost 30% of the minimum sepa- 
ration claimed by Coles et al. (2004). This emphasizes 
the importance of incorporating lensing effects into the 
analysis of scintillation correlations which could further 
refine the estimate of inclination obtained by Coles et al. 
(2004). 

Light bending also causes the radio beam of pulsar 
A to pass through the regions in the pulsar B mag- 
netosphere which would otherwise not be sampled by 
the beam. Since the magnetospheric absorption depends 
strongly on the distance from pulsar B (Rafikov & Gol- 
dreich 2005), one expects gravitational lensing to have an 
important effect on the eclipse profile of pulsar A. Varia- 
tions in the pattern of magnetospheric absorption would 
then be expected on a timescale of apsidal precession be- 
cause of the change in the minimum projected distance 
between the two pulsars at conjunction (see H2.4(l . 

Pulsar A can produce gravitational lensing of pulsar 
B in exactly the same way as B does it with pulsar A. 
The Einstein radius of pulsar A is slightly larger than 
that of pulsar B, 2640 km instead of 2550 km, and the 
measurements would probably be cleaner in this case (al- 
though the timing accuracy of B is not as good as that 
of A) since the radio beam of pulsar B would not pass 
through the magnetosphere of pulsar A (with light cylin- 
der radius Rl,a ~ 1000 km). Unfortunately, pulsar B is 
currently very weak at the orbital phases when it is be- 
hind pulsar A (Ramachandran et al. 2004) which makes it 
hard to detect lensing in this orbital configuration. How- 
ever, the system's geometry will change due to the effects 
of general relativity and pulsar B might become active 
and bright at its superior conjunction (Jenet & Ransom 
2004). This may give us additional handle on the grav- 
itational lensing in J0737 and permit the determination 
of the pulsar B spin axis orientation (see ij2.1|l . 
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effect of the plasma even less important. 
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